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Abstract 

We propose a formal extension of thermodynamics and kinetic theories to a larger 
class of entropy functionals. Kinetic equations associated to Boltzmann, Fermi, Bose and 
Tsallis entropies are recovered as a special case. This formalism first provides a unifying 
description of classical and quantum kinetic theories. On the other hand, a generalized 
thermodynamical framework is justified to describe complex systems exhibiting anoma- 
lous diffusion. Finally, a notion of generalized thermodynamics emerges in the context 
of the the violent relaxation of collisionless stellar systems and two-dimensional vortices 
due to the existence of Casimir invariants and incomplete relaxation. A thermodynamical 
analogy can also be developed to analyze the nonlinear dynamical stability of stationary 
solutions of the Vlasov and 2D Euler-Poisson systems. On general grounds, we suggest 
that generalized entropies arise due to the existence of "hidden constraints" that modify 
the form of entropy that we would naively expect. Generalized kinetic equations are there- 
fore "effective" equations that are introduced heuristically to describe complex systems. 

1 Introduction 

Standard kinetic equations satisfy two fundamental properties linked to the first and second 
principles of thermodynamics: the conservation of energy (and mass) and the increase of entropy 
(H-theorem). These properties are shared in particular by the Boltzmann and by the Landau 
equations which are at the basis of the kinetic theory of dilute gases and neutral plasmas pQ . 
When the system is in contact with a thermostat, instead of being isolated, the conservation of 
energy and the increase of entropy (microcanonical description) are replaced by the decrease 
of free energy F = E — TS at fixed temperature (canonical description). This is the proper 
description of Brownian motion that is analyzed in terms of stochastic processes and Fokker- 
Planck equations (Kramers, Smoluchowski,...) |2j. In the standard framework, the functional 
increasing monotonically with time is the Boltzmann entropy Ss[f] — ~ J f In fd 3 rd 3 v or the 
Boltzmann free energy JbU] = SbII] ~ /3E[f}. 

In a recent paper |3] , we have proposed to develop a generalized thermodynamical formalism 
for a larger class of functionals that we called generalized entropies. They can be written 
S[f] — — J C(f)d 3 rd 3 v where C(f) is a convex function, i.e. C"(f) > 0. Boltzmann, Fermi, 
Bose, and Tsallis entropies are particular functionals of the above form. On general grounds, 
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what we mean by generalized thermodynamics is the extension of the usual variational principle 
of classical thermodynamics (maximization of Boltzmann entropy at fixed mass M and 

energy E) to a larger class of functionals. This variational principle arises in many problems 
of physics (or biology, economy,...) for various reasons that do not have necessarily a direct 
relation to thermodynamics. In many cases, it is relevant to develop a thermodynamical analogy 
and to use the same vocabulary as in standard thermodynamics. This allows one to transpose 
directly the standard methods developed in ordinary thermodynamics to a new context. 

In [3] we have proposed a generalized class of Fokker-Planck equations associated to this 
generalized thermodynamical framework. This formalism is interesting to develop because it 
leads to a unified description of known kinetic theories (classical and quantum) and it also 
generates new types of kinetic equations. In that respect, it can be of interest in applied math- 
ematics and theoretical physics. This formalism can also have important physical applications. 
For example, these generalized Fokker-Planck equations can account for a process of anomalous 
diffusion in complex media. They arise naturally from ordinary Fokker-Planck equations by 
assuming that the diffusion coefficient is a function of the density. Generalized thermodynam- 
ics can also be relevant for the violent relaxation of stellar systems and two-dimensional (2D) 
vortices. In that context, generalized entropies (also called if -functions) emerge due to the 
existence of fine-grained constraints (Casimir invariants) that modify the form of entropy that 
we would naively expect. Generalized Fokker-Planck equations can provide a simple small- 
scale parametrization of turbulence (mixing) in stellar dynamics and 2D hydrodynamics. They 
can also serve as powerful numerical algorithms to compute arbitrary nonlinearly dynamically 
stable solutions of the 2D Euler-Poisson or Vlasov-Poisson systems. Indeed, the condition of 
nonlinear dynamical stability can be put in a form analogous to a condition of generalized 
thermodynamical stability jUEllHllZj- This is a striking illustration of the thermodynamical 
analogy mentioned above. Since the notion of generalized thermodynamics can have different 
interpretations, it is relevant to work at a general level and develop a formalism without ex- 
plicit reference to a precise context. Then, a justification of this generalization and a physical 
interpretation of the results must be given in each case. 

In an early work [H], we observed that classical and quantum Fokker-Planck equations can 
be obtained from a phenomenological Maximum Entropy Production Principle (MEPP) 
by maximizing the rate of entropy (resp. free energy) production at fixed mass and energy 
(resp. temperature). This variational approach, closely related to the linear thermodynamics 
of Onsager, is the most natural extension of the equilibrium thermodynamical principle in which 
we maximize entropy (resp. free energy) at fixed mass and energy (resp. temperature). In [3], 
we generalized this principle to a larger class of entropy functionals and obtained generalized 
Fokker-Planck equations. We also "guessed" a form of generalized Landau equation consistent 
to the generalized Kramers equation obtained with the MEPP. Recently, we found that a 
similar generalization of kinetic theory was attempted by Kaniadakis ^0] from a different point 
of view. His approach consists in generalizing the assumptions that are made at the start to 
derive the Boltzmann and the Fokker-Planck equations. This amounts to modifying the form 
of the transition probabilities that arise in the dynamical process. Such generalized transition 
probabilities are relevant for quantum particles (fermions and bosons) and possibly, also, in the 
physics of complex media. The MEPP approach, on the other hand, is purely thermodynamical 
and exploits at best the first and second principles of thermodynamics (possibly extended to 
generalized functionals) in a viewpoint reminiscent of Jaynes' ideas. 

In this paper, we show that the two approaches lead to equivalent kinetic equations. In Sec. 
121 we use Kaniadakis approach to derive the generalized Landau equation from the generalized 
Boltzmann equation in a weak deflexion approximation. In Sec. OJ we show that the generalized 
Landau equation can also be obtained by coarse-graining the Vlasov equation in the context 
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of the violent relaxation of stellar systems. In Sec. 0] we establish the main properties of 
the generalized Landau equation (conservation laws, generalized H-theorem,...). In Sec. EJ 
we derive the generalized Kramers equation from the generalized Landau equation in a test 
particle approach and a thermal bath approximation. The generalized Smoluchowski equation 
is in turn derived from the generalized Kramers equation in a hydrodynamical limit. In Sees. El 
and we determine explicit expressions of the diffusion coefficient for Boltzmann, Fermi and 
Tsallis distributions. In Sec. |Hl we use the generalized Kramers equation to derive generalized 
truncated distribution functions accounting for an escape of particles above a limit energy. 

2 Generalized kinetic equations 

2.1 The generalized Boltzmann equation 

The ordinary Boltzmann equation can be written in the form 



where dfl is the element of solid angle and w(v, Vi; v', v^) is the density probability of a collision 
transforming the velocities v,vi in v',v[ or the converse (the abbreviations jf(v'), /(vj_), /(v), 
/(vi) stand for for /(r, v', t), f(r,v[,t), f(r,v,t), /(r, v', i)). The material derivative is d/dt = 
d/dt + vd/dr + Fd/dv where F = — V$ is a mean-field force acting on the particles. The 
Boltzmann equation conserves the mass 



(1) 




(2) 




the energy 



(3) 




the angular momentum 



(4) 




and the impulse 



(5) 




where p(r, t) is the spatial density. In addition, the Boltzmann entropy 



(6) 




satisfies a if -theorem, i.e. S > with S = if, and only if, the distribution /(r, v) is the 
Boltzmann distribution 



(7) 



feq(r, v) = Ae 



where e' = y 



+ $ +Q • (rxv) + U • v is the energy of a particle by unit of mass. 



3 



Recently, Kaniadakis [TU] has proposed the generalization 

(8) | = | d 3 vi dn w(v, Vl ; v', vi)|a(/')6(/)a(/i)6(/i) - a{f)b{f)<fiWi) 

where the functions a(f) and b(f) are somewhat arbitrary (we have noted / = /(v), j\ = f(yi), 
f = /(v') and /j = /(v^)). This generalization encompasses the case of quantum particles 
(fermions and bosons) with exclusion or inclusion principles. This generalization could also be 
relevant in the case of complex systems for which the transition probabilities are not simply the 
product of distribution functions. This can happen when we are not in the strict conditions of 
validity of the ordinary Boltzmann equation. It is also possible that Eq. (jHJ) is just an effective 
kinetic equation accounting for "hidden constraints" in complex media. 

In the following, we shall consider the situation in which the potential of interaction between 
particles is Coulombian (or Newtonian). In that case, each encounter provokes a weak deflexion 
of the particles trajectory and it is of order to consider the weak deflexion limit of the Boltzmann 
equation. Classically, this leads to the so-called Landau equation which forms the basis of the 
kinetic theory of neutral plasmas [Tj and stellar systems [TT]. Our aim is this section is to 
derive a generalized Landau equation from the generalized Boltzmann equation proposed by 
Kaniadakis. This generalization is essentially formal. In the following, we follow the classical 
derivation of the Landau equation reported in the monograph of Balescu pQ. Therefore, we shall 
omit the calculations that are identical to the classical case and refer to |X| for more details. 

2.2 The weak deflexion approximation 

First, it is convenient to write the velocities of the particles before and after the collision as 

(9) v' = v + A, 



(10) v£ = vi - A, 

where A is the velocity deviation. We can now express the probability of a collision in terms 
of new variables as P: 

A A 

(11) w(v, vi; v', vi) — > w(v + — , vi — — ; A). 
The generalized Boltzmann equation can thus be rewritten 



dt 



J d 3 vi dn w(v + ^ , vi - y ; A) ja[/(v + A)]6[/(v)]o[/(vi - A)]6[/(vi)] 
12) -aifWmv + AMHvxMf^-A)]). 



In the weak deflexion limit |A| <C |v|, |vi|, we can expand the r.h.s. of Eq. (|12j) in Taylor 
series. Using 

. . .A A AN , . 1 A „ (dw(v, vi; A) dw(v, Vi; A) 1 

(13) -(v + T , v, - —; A) = „(v, v i; A) + ^A^j QyP -tr^\ + 

df(v) 1 d 2 f(v) 

(14) /(v + A) = /(v) + A^H^ + -A^A" J{ ' 1 



dvv 2 dv^dv v 
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we get 



| = \ J d 3 Vl dOA'A^v, Vl ; A) [aKMi - + abifra" - ^f|f| 



<9 2 / dfdfi dfdf 

+ai6i (6a' — 6'a) — — 2 (661 a'a'i — aaib'b'-,)— — — h ai6i (6a" — ab")— — — — 

dv^dv" dv^ dv\ ;>,>/<;),.'' 



dw 



(15) 



dw 



9f 9fi 
a\b\(ba! — 6'a) — ab(bia' x — a\b'i 



dv l 



dv\ 



0161(60' - b'a) — - 06(6101 - a i & 'i)^ 



where a = a[/(v)], ai = a[/(vi)], a' = a'[/(v)], a' x = o'[/(vi)] etc... Integrating the last term 
by parts, the foregoing expression simplifies in 



d£ = 1 
dt 2 

(16) 



d 3 Vl d^A^A^ w(v, vi; A 



<9 2 / dfdfi 
0161(60' — 6'a)— — (6ia^ — 6' 1 ai)(a6' + a'6) 



dv^dv u 



+aibi(ba" — a6" 



9f ^ <9t> " 



+ 



<9«; 



df 1 ' / 1 

ai6 x (6a' - h ' a )^ ~ ob( h i a 'i - aib '^g^ 



Equation fjlCjf) can be written more compactly as 

(17) I = ^ rf 3 Vi ^| ai6l(W _ ^ _ a6(M ; _ 



where we have defined 

(18) K>* U = ^J dQ w(v, vi, A)A M A I/ . 
We now define the functions g and o by 

(19) 9(f) = a(f)b(f), h(f) = b(f)a'(f) - b'(f)a(f). 
Then, Eq. (117(1 can be rewritten 



(20) 



dt 



rfV K^{g{h)h{f)^-g{f)h{h) d J± 



This will be called the generalized Landau equation. The tensor K^ v can be calculated explicitly 
in the linear trajectory approximation pQ. The result can be expressed as 



(21) 



A 



( 



5^ 



u \ 



where u = Vi — v is the relative velocity and A is a constant (in the plasma case A = 
(e 4 /8iim 3 el)\ia(L Deb y e /L min ) and in the gravitational case A = 2nG 2 m\a.(L max / L min )). 

By developing a kinetic theory of 2D point vortices |T21 EE! we have derived a kinetic 
equation of the form 



(22) 



dt w 8 <9r^ 



/ dWK .v)(p^-P^ 
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(23) K'"' 



where 

and £ = Yi — r, v = (V)(r 1; t) — (V)(r, t). This equation conserves all the constraints of the 
point vortex model and increases the Boltzmann entropy (if-theorem). It is reminiscent of the 
Landau equation but it differs from the Landau equation due to the 5-function which takes 
into account the conservation of energy E = \ J utpcPr where ip(r,t) is the stream-function. 
In addition, its physical interpretation and derivation is completely different from that of the 
Landau equation. We can heuristically propose a formal extension of this equation to the more 
general form 

dP Ni 2 d 

(24) ^r + (V)VP - 



dt 8 dr^ 



a. 



2.3 Generalized entropy 

We shall say that a kinetic equation possesses a generalized microcanonical thermodynamical 
structure if it conserves mass and energy and increases continuously a functional of the form 

(25) S = - J C(f)dhd\, 

where C(f) is a convex function. By analogy with ordinary thermodynamics, the functional 
(|23j) will be called a generalized entropy even if it does not correspond to a true entropy in the 
strict sense (in that case, we speak of a thermodynamical analogy). The equilibrium distribution 
function f eq {r, v), reached by the kinetic equation for t — > +oo, maximizes S[f] at fixed E and 
M. Introducing Lagrange multipliers and writing the variational principle in the form 

(26) 5S - (35E - a5M + (3S15L + (3V5P = 0, 

where we have also accounted for the conservation of angular momentum and impulse, we find 
that the equilibrium distribution function is given by 

(27) C(f eq ) = -p(j + $) +/3ft(rxv)+/3Uv 

As we shall see, the generalized Landau equation ()20j) possesses a microcanonical thermo- 
dynamical structure. The generalized entropy can be expressed in terms of the functions a(f) 
and b(f) as 

(28) S = - J C(f)d 3 rd\, with C\f) = In ^ 

From Eqs. (|28j) and (|19jh we obtain the relation 

(29) h(f)=g(f)C"(f), 
which leads to the identities 

(30) a(f) = vW)e^' (/) , &(/) = vW)e^ C ' (/) - 

If we know a and b, we can obtain g and h from Eq. (fTT?|) and C from Eq. ()28j) . Alternatively, 
if we know h and g, we obtain C"(f) from Eq. (p§|) and deduce a and b from Eq. (|5U|) (up to a 
multiplicative factor). If we only specify C(f), we cannot obtain a and b individually but only 
the ratio a/b. In the next section, we shall consider simplified forms of the generalized Landau 
equation where everything is determined by the specification of the generalized entropy S[f], 
or equivalently by the function C(f). 
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2.4 Simplified forms of the generalized Landau equation 

We shall first impose that 

(31) 9(f) = f, h(f) = fC"(f), 
where C is a convex function. In that case, Eq. (J2(Jj) takes the form 

(32) |^/^^4 C " (/, ^- C " (/l) i}- 

This generalized Landau equation was first written in j3J. From Eqs. (|3Uj) and ()31|) we find 
that a(/) and b(f) are related to C(/) by 

(33) a(f) = y/feW, b(f) = Jfe'^. 

It is interesting to consider particular cases of the generalized Landau equation (|32|). For the 
Boltzmann entropy C(f) = / In/, we recover the ordinary Landau equation 

(34) * '/^W,* 



(it y \ dv v dv\ 

which is the weak deflexion limit of the ordinary Boltzmann equation corresponding to 
(35) a(f) = /, b(f) = 1, 

in Eq. (jSJ). For the Tsallis entropy C(f) = — ^-(Z 9 — /), we obtain the g-Landau equation 

which is associated to the g-Boltzmann equation (jHJ) with 

(37) o(/) = vZ/e^T^" 1 - 1 ), &(/) = v^e-^i)^ 9 "- 1 ). 
Instead of Eq. (jSH), we can impose the relations 

(38) h(f) = 1, g(f) = 1/C"(f). 
In that case, Eq. (j2U)l reduces to 

(39) d f = °[ d s VlK ^ 1 9/ 1 9/1 



eft I 9u" C"(f)dv1 

This alternative form was also written in [Hj. From Eqs. (jHDj) and (|38jl . we find that a(/) and 
b(f) are related to C(/) by 



(40) a(f) = -_^ e ^'(/) 6 (/) = _^ =e -^'(/)_ 

For the entropy C(f) = / In / + (r) - fif) \n(r] - [if), we get 



For /i — 1, the foregoing equation describes the case of fermions accounting for the Pauli 
exclusion principle. For // = 1, it describes the case of bosons. For intermediate values of /i 
it describes a intermediate quantum statistics interpolating between the Bose and Fermi ones. 
The generalized Landau equation ()41|) is the weak deflexion limit of the generalized Boltzmann 
equation (jSJ) with 

(42) a(f) = f, b(f) = Vo - nf. 

Although this extension of kinetic theory is interesting on a formal point of view, we do 
not claim that it is relevant to plasma physics and stellar dynamics (except in the quantum 
case). Indeed, in neutral plasmas where the interaction is short-ranged due to Debye shielding 
the diffusion is normal and the ordinary Landau equation is rigorously valid. On the other 
hand, in collisional stellar systems, the diffusion is only slightly anomalous due to logarithmic 
divergences p%j . Therefore, the ordinary Landau equation remains marginally valid when 
logarithmic divergences are properly regularized [T^]. There can be corrections to the ordinary 
Landau equation due to memory effects and spatial derealization (similar effects arise in 
the kinetic theory of point vortices [121 )• However, this does not apparently justify a rigorous 
notion of generalized thermodynamics (even if deviations to the Maxwellian distribution may 
be expected for intermediate times). For large times, ordinary thermodynamics (based on the 
Boltzmann entropy) is rigorously justified for collisional stellar systems (and point vortices) 
in a suitable thermodynamic limit although these systems are non-extensive and non-additive 
(see for a more complete discussion). Therefore, there does not seem to be any justification 
of a generalized thermodynamics for Hamiltonian systems of point particles in the infinite 
time limit (collisional relaxation). At the present time, it is not clear to which systems the 
generalized kinetic theory developed previously could apply (with the exception of quantum 
particles). This is an open problem left for future investigations. Our guess is that generalized 
kinetic equations can serve as effective equations in the case of complicated systems. In the 
following section, we show that a notion of generalized thermodynamics also emerges in the 
context of the violent relaxation of collisionless stellar systems (and other Hamiltonian systems 
with long-range interactions) for a completely different reason. Generalized kinetic equations 
can find physical applications in that context. 



3 Violent relaxation of collisionless stellar systems 
3.1 The Vlasov-Poisson system 

For most stellar systems, including the important class of elliptical galaxies, the encounters 
between stars are completely negligible [TT1 [T7] and the galaxy dynamics is described by the 
self-consistent Vlasov-Poisson system 



(44) A$ = 4ttG J fd 3 v. 

Here, /(r, v,t) denotes the distribution function (defined such that fd 3 rd 3 v gives the total 
mass of stars with position r and velocity v at time t), F(r, t) = — V$ is the gravitational 
force (by unit of mass) experienced by a star and $(r,t) is the gravitational potential related 
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to the star density p(r, t) — J fd 3 v by the Newton-Poisson equation (|44j) . The Vlasov equation 
PHj) simply states that, in the absence of encounters, the distribution function / is conserved 
by the flow in phase space. This can be written df/dt = where d/dt = d/dt + UeV6 is the 
material derivative and XJq = (v, F) is a generalized velocity field in the 6-dimensional phase 
space (r, v) (by definition, V6 = (d/dr, d/d~v) is the generalized nabla operator). Since the flow 
is incompressible, i.e. VeU 6 = 0, the hypervolume of a "fluid" particle is conserved. Since, in 
addition, a fluid particle conserves the distribution function, this implies that the total mass (or 
hypervolume) of all phase elements with phase density between / and f + Sf is conserved. This 
is equivalent to the conservation of the Casimir integrals Ih = J h(f)d 3 rd 3 v for any continuous 
function h(f). This is also equivalent to the conservation of the moments 

(45) M n = J f n d 3 rd 3 v, 

which include in particular the total mass M = j fd 3 vd 3 \. It is also straightforward to check 
that the Vlasov-Poisson system conserves the total energy E, the angular momentum L and 
the impulse P. 



3.2 The metaequilibrium state 

The Vlasov-Poisson system develops very complex filaments as a result of a mixing process 
in phase space. If we introduce a coarse-graining procedure, the coarse-grained distribution 
function /(r,v,t) will reach a metaequilibrium state f(r, v) on a very short timescale, of the 
order of the dynamical time. This process is known as "phase mixing" and "violent relaxation" 
|1 lj . Lynden-Bell ^Bj has tried to describe this metaequilibrium state in terms of statistical 
mechanics. If p(r, v, 77) denotes the density probability of finding the value / = 77 of distribution 
function in (r, v), then the mixing entropy is given by 

(46) S[p} = - J p In p d 3 rd 3 vdr]. 

It can be obtained by a standard combinatorial analysis [IE] or by using the concept of Young 
measures and large deviations [19 1. Assuming ergodicity (which may not be realized in practice, 
see below) the statistical equilibrium state is obtained by maximizing S[p] while conserving mass 
M, energy E and all the Casimirs (or moments M n ). The optimal equilibrium state can be 
written [TBI [7j 

(47) p ( r)V;r/ ) = I x ( r? )e-(/ 3 ^, 

where e = ^- + $ is the energy of a star by unit of mass, x(v) = exp(— J2 n >i a nV n ) accounts 
for the conservation of the fragile moments M n> \ = j pi] n dr]d 3 rd 3 v and a, (3 are Lagrange 
multipliers for M and E (robust integrals). The partition function Z is determined by the local 
normalization f pdr\ = 1 and the equilibrium coarse-grained field is given by / = J pr\dr\. This 
can be written [2] 

(48) 7=~^f = F((3e + a)=J(e). 

We emphasize that Eq. (J48|) is obtained after two successive coarse-grainings. Intrinsically, 
a galaxy is a collection of iV point-like stars and the exact distribution function f exact = 
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X^=i m «<K r — r i)^( v — v «) is a sum °f (^-functions. As is customary in statistical mechan- 
ics, we consider a smooth distribution function f(r,v,t) which is the statistical average of 
f exact, i-e. / = (f exact)- In the collisionless regime valid for times t <C t coU ~ 

In N^dyn (tcoll IS the 

timescale of collisional relaxation and td yn the dynamical time), this function /(r, v,t) satisfies 
the Vlasov equation ()43|1. However, the Vlasov-Poisson system develops a mixing process and 
it is relevant to introduce another coarse-graining (at a much larger scale) and define f(r,v,t) 
as the local average of /(r, v, t) on a phase space cell of volume e^e 3 ,. The quantity appearing in 
Eq. (|4"8j) is the most probable form of / at statistical equilibrium (in the sense of Lynden-Bell), 
assuming ergodicity (i.e. efficient mixing). 

From Eq. (jSj) , it is easy to show (see the equivalent proof in [20]) that 

(49) f(e) = -0f 2 , h = J P(V ~ 7) 2 dv > 0, 

where f 2 is the centered local variance of the distribution p(r,v,rj). Therefore, / = /(e) is 
a decreasing function of the stellar energy (assuming (3 > 0). Since /(e) is monotonic, the 
coarse-grained distribution ([4*%)) extremizes a "generalized entropy" |7J |3] 

(50) S[f] = -J C(/)rf 3 rrf 3 v, 

at fixed mass M and energy E, where C(f) is a convex function, i.e. C" > 0. Indeed, 
introducing Lagrange multipliers and writing the variational principle in the form 

(51) 5S - f35E - a5M = 0, 
we find that 

(52) C'(f) = -Pe-a. 

Since C is a monotonically increasing function of /, we can inverse this relation to obtain 

(53) f = F(/3e + a)=f{e), 

where F(x) = (C") _1 (— x). Equation (jHSj) can be compared to Eq. pSjl. From the identity 

(54) T{e) = -P/C"(f), 

resulting from Eq. ()52|). /(e) is a monotonically decreasing function of energy if (3 > 0. There- 
fore, for any Gibbs state of the form (|47|). there exists a generalized entropy of the form (|50|) 
that / extremizes (at fixed E, M). It can be shown furthermore that / maximizes this func- 
tional. We note that C(f), hence the generalized entropy (|50|1. is a non-universal function 
which depends on the initial conditions. In general, it is not the ordinary Boltzmann entropy 
Sb[J] — — J f hi fd 3 rd 3 v due to fine-grained constraints (Casimirs) that modify the form of 
entropy that we would naively expect. We emphasize that maximizing the multi- levels Boltz- 
mann entropy S[p] at fixed mass M, energy E and with an infinite number of fine-grained 
constraints M n (Casimirs) gives the same result for the coarse-grained field / as maximizing a 
certain generalized entropy S[f] (non-universal) while conserving only mass M and energy E 
(robust constraints). The existence of "hidden constraints" (here the Casimir invariants that 
are not accessible at the coarse-grained scale) is the physical reason for the occurence of "gen- 
eralized entropy" functionals in a problem. We can either work with the Boltzmann entropy 
and take into account all the constraints imposed by the dynamics, or keep only the constraints 
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that are the most directly accessible to the observations and change the form of entropy. This 
clearly leads to an indetermination which appears in the parameter q of Tsallis or more gener- 
ally in our function C(f). We expect this idea of "hidden constraints" to be very general and 
of fundamental importance. 

The above statistical approach rests on the assumption that the evolution is ergodic. In 
reality, this is not the case. It has been understood since the beginning ^H] that violent 
relaxation is incomplete so that the Boltzmann entropy (pfUj) is not maximized in the whole 
available phase space (this is independant on the fact that it has no maximum!). However, the 
metaequilibrium state reached by the system as a result of incomplete violent relaxation is always 
a nonlinearly dynamically stable solution of the Vlasov-Poisson system (on the coarse-grained 
scale). If / = /(e), which is a particular case of the Jeans theorem [H], it maximizes a functional 
of the form ()5(J|) at fixed mass and energy. In this dynamical context, S[f] is called a if -function 
j^]. This functional depends on the initial conditions (for the same reasons as before) and also on 
the strength of mixing (if the mixing is complete, C(f) can be predicted by the statistical theory 
of Lynden-Bell). Boltzmann and Tsallis functionals are particular H-functions corresponding 
to isothermal and polytropic distribution functions. They are not good models of incomplete 
relaxation for elliptical galaxies [7] . A better model is a composite model that is isothermal in 
the core and polytropic in the halo. Since the variational principle determining the nonlinear 
dynamical stability of a collisionless stellar system (maximization of a H-function at fixed mass 
and energy) is similar to the usual thermodynamical variational principle (maximization of the 
Boltzmann entropy at fixed mass and energy) we can use a thermodynamical analogy to analyze 
the dynamical stability of collisionless stellar systems E] . In this analogy, the H-function 
can be called a "generalized entropy". We believe that this thermodynamical analogy is the 
correct interpretation of the notion of "generalized thermodynamics" introduced by Tsallis in 
the context of stellar systems (and 2D turbulence). We emphasize that the maximization of a 
if -function (e.g., Tsallis entropy) at fixed mass and energy is a condition of nonlinear dynamical 
stability, not a condition of thermodynamical stability. 

3.3 Heuristic approach of violent relaxation 

Violent relaxation is a very complicated concept because of the presence of fine-grained con- 
straints (Casimirs). The Casimirs differ from robust constraints such as mass and energy 
because they are altered by the coarse-graining procedure since f n ^ f for 1. Therefore, 
they can be determined only from the initial conditions (which are not mixed) at t = 0, say. 
Unfortunately, in practice, we do not know the initial conditions (e.g., the initial state that 
gave rise to an elliptical galaxy) so that we do not know the Casimirs. Only mass and energy 
(robust constraints) can be determined at all times t > since M = M and E ~ E . We thus 
have to deal with a limited amount of information on the system. Therefore, we cannot make 
predictions because we do not know all the constraints on the system. If we want to make some 
predictions, we have two possibilities: (i) the first possibility is to "guess" the initial conditions 
(consistent with the information that we have on the system) and determine the corresponding 
equilibrium state. We can then study how the equilibrium state depends on the initial con- 
ditions (for given E and M). (ii) the second possibility is to "guess" the generalized entropy 
that is maximized by the system at equilibrium. We can then study how the equilibrium state 
depends on the form of C(f) and whether generalized entropies presenting the same properties 
can be regrouped in "classes of equivalence" 0. 

Each approach has its own advantages and drawbacks. The first approach is more compli- 
cated because, for realistic initial conditions, we have to account for an infinity of constraints 
(Casimirs) in addition to mass and energy. This clearly leads to practical difficulties. In ad- 



11 



dition, it is not clear whether all these constraints are physically relevant because they can 
be altered by non-ideal effects as discussed in [7J Ej. Finally, this description assumes that 
the mixing is complete which is not the case in practice. We shall therefore prefer the second 
approach where we have only two robust constraints M and E. The other constraints are taken 
into account implicitly in the form of the generalized entropy that we consider. This is a much 
more convenient approach. In addition, this approach is directly connected to the problem of 
dynamical stability that we discussed previously. It has therefore a lot of attractive advantages. 



3.4 The quasilinear theory 

Basically, a collisionless stellar system is described in a self-consistent mean field approxima- 
tion by the Vlasov-Poisson system (JJ3J) (JHJ) . In principle, these coupled equations determine 
completely the evolution of the distribution function /(r, v, t). However, as discussed in Sec. 
13.21 we are not interested in practice by the finely striated structure of the flow in phase space 
but only by its macroscopic, i.e. smoothed-out, structure. Indeed, the observations and the nu- 
merical simulations are always realized with a finite resolution. Moreover, the "coarse-grained" 
distribution function /(r, v, t) is likely to converge towards an equilibrium state contrary to the 
exact distribution / which develops smaller and smaller scales. 

If we decompose the distribution function and the gravitational potential in a mean and 
fluctuating part (/ = / + /,$ = $ + $) and take the local average of the Vlasov equation ^ 
we readily obtain an equation of the form 

d l d 7 t^J dJj 
(55) + v^- + F- 



dt dr dv dv 

for the "coarse-grained" distribution function with a diffusion current J/ = /F related to 
the correlations of the "fine-grained" fluctuations. Any systematic calculation of the diffusion 
current starting from the Vlasov equation (|43|) must necessarily introduce an evolution equation 
for /. This equation is simply obtained by subtracting Eq. (}55|) from Eq. (J43|) . This yields 



(56) | + 4 +P w = _*«_ f w +f v 

dt dr dv dv dv dv 

To go further, we need to implement some approximations. In the sequel, we shall develop 
a quasilinear theory (see also [5111221 El)- This will provide a precise theoretical framework to 
analyze the process of "collisionless relaxation" in stellar systems. The essence of the quasilinear 
theory is to assume that the fluctuations are weak and neglect the nonlinear terms in Eq. (|56j) 
altogether. In that case, Eqs. (J55j) and reduce to the coupled system 

dl - d — 



where L = v-j^ + F^ is the advection operator in phase space. Physically, these equations 
describe the coupling between a sub dynamics (here the small scale fluctuations /) and a macro- 
dynamics (described by the coarse-grained distribution function /). Due to the strong simpli- 
fications implied by the neglect of nonlinear terms in Eq. (|56jh the quasilinear theory only 
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describes the late quiescent stages of the violent relaxation process, when the fluctuations have 
weaken (gentle relaxation). Although this is essentially an asymptotic theory, it is of impor- 
tance to develop this theory in detail since it provides an explicit expression for the effective 
"collision operator" which appears on a coarse-grained scale. 
Introducing the Greenian 

(59) G(h,ti) =exp|-^ dtL(t) 
we can immediately write down a formal solution of Eq. (J58j) , namely 

(60) /(r, v, t) = Q(t, 0)/(r, v, 0) - J dsG(t, t - s)F(r, t - s)S-(r, v,t-s). 

Although very compact, this formal expression is in fact extremely complicated. Indeed, all 
the difficulty is encapsulated in the Greenian Q(t,t — s) which supposes that we can solve the 
smoothed-out Lagrangian flow 

/ % dr dw — 

(61) M =¥ ' 

between t and t — s. In practice, this is impossible and we will have to make some approxima- 
tions. 

The objective now is to substitute the formal result (jHUj) back into Eq. (JH7|) and make 
a closure approximation in order to obtain a self-consistant equation for f(r, v, £). If the 
fluctuating force F were external to the system, we would simply obtain a diffusion equation 

(62) S. +L J-JL(d^K 
with a diffusion coefficient given by a Kubo formula 

(63) = [ dsF»(r,t)F v (r,t-s). 

Jo 

However, in the case of the Vlasov-Poisson system, the gravitational force is produced by the 
distribution of matter itself and this coupling will give rise to a friction term in addition to the 
pure diffusion. Indeed, we have 

(64) F(r,t) = j F(r' ^r)f(r',v',t)d 3 r'd 3 v', 
where 

( 65 ) F (r - = G ^Ts> 



r — r 



represents the force (by unit of mass) created by a star in r' on a star in r (Newton's law). 
Therefore, considering Eqs. (|5U|) and (jEU), we see that the fluctuations of the distribution 
function f(r,v,t) are given by an iterative process: f(t) depends on F(t — s) which itself 
depends on f(t — s) etc... We shall solve this problem perturbatively in an expansion in powers 
of the gravitational constant G. This is the equivalent of the "weak coupling approximation" 
in plasma physics. 
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According to Eq. (|64|). we have 
(66) /(r, v, t)F»(r, t) = J d 3 r'd 3 v'F(r' r)/(r', v', £)/(r, v, £). 

On the other hand, according to Eqs. (JoUj) and (p>4"j) . we have 



/(r, v, t) = 0(t, 0)/(r, v, 0) - J ds J d d r"d 6 v"g(t, t - s) 

(67) xF"(r" -> r)/(r", v",t - s)^(r, v,t - s), 

and a similar expression for f(r',v',t). Substituting the foregoing expansion in Eq. ()66|) . we 
find that the term of order G 2 is 

f{T,v,t)F»{T,t) = - J^ds J d 3 r"d 3 v" J d 3 r'd 3 v'F^(r' -> r) 
x 0)£(t, t - s)F»(r" - ^/(r^O^r^vV-s)^^, v, t - s) 

= 97 

(68) 0)fiT(t, t - S )F"(r" - r')/(r, v, 0)/(r", v", t - s)-^(r', v', t - s) 
Using 

(69) g , (t,o) = g%t-s)g , (t-s,o), 

and 

(70) /(r', v', t - a) = f?'(f - a, 0)/(r', v', 0) + 0(G), 
we obtain after some rearrangements 



+ Lf = -^ J^ds j dVdVcfYW'F'V -> r)0'(t, t - s)0(*, t - s) 
x |^(r" - r)/(r', v', t - s)/(r», v", t - s)J^(r, v, t - s) 
+J F V - r')/(r, v, t - s)f(r", v", t - s)^^, v', t - s) 



(71) 

In this expression, the Greenian Q refers to the fluid particle r(£), v(t) and the Greenian Q' to 
the fluid particle r'(i), v'(i). To close the system, it remains for one to evaluate the correlation 

function f(r,v,t)f(r',v',t). We shall assume that the mixing in phase space is sufficiently 
efficient so that the scale of the kinematic correlations is small with respect to the coarse- 
graining mesh size. In that case, 



(72) /(r, v, t)f(r>, v', t) = e 3 e 3 v 5(r - r')5(v - v')/ 2 (r, v, t), 

where e r and e v are the resolution scales in position and velocity respectively and 



(73) h = P = (/ - /) 2 = P ~ f- 

is the local variance of the fine-grained fluctuations. 
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3.5 The closure approximation 

We are now led to a closure problem. Indeed, in order to obtain a self-consistent kinetic 
equation for /, we need to determine the variance f%. If the initial condition in phase space 
consists of patches of uniform distribution function f = rj surrounded by vacuum / = (two- 
levels approximation), then f 2 = r] xf = rj f and, therefore, / 2 = f(r) — f). This leads to a 
generalized Landau equation of the form (|%Tj) as for fermions 122 CZ] • For more complicated 
initial conditions (multi- levels case), the strategy would be to write down a kinetic equation for 
p(r, v, rj, t), the density probability of finding the value / = rj in (r, v) at time t. This extension 
can be performed along the lines sketched in J22J and the resulting equation for p(r, v, rj, t) can 
be closed. However, this approach leads to a system of N coupled equations (one for each level 
77) which is not convenient to solve when N ^> 1. We could alternatively obtain a hierarchy 
of equations for the moments /, f 2 ,...,f n ,... but we would then encounter a closure problem. 
The equations obtained by this method are complicated because they take into account the 
conservation of all the Casimirs. 

In this paper, we propose a closure approximation that leads to a simpler kinetic equation. 
While not being exact, this equation preserves the robust features of the process of violent 
relaxation and is amenable to an easier numerical implementation. Its main interest is to go 
beyond the two-levels approximation while leaving the problem tractable. The idea is to observe 
that Eqs. (jjHj) and (J53)l lead to the important relation 

(74) f 2 



C"(f) 

This relation is valid at equilibrium but we propose to use it as a closure approximation in 
Eq. (|72|). This is expected to be a reasonable approximation if we are close to equilibrium, 
which is in fact dictated by the quasi-linear approximation. Of course, this procedure assumes 
that we know the function C(f) in advance. This is the case if we have already determined 
the equilibrium state and we want to describe the dynamics close to equilibrium. This is also 
the case in the heuristic approach of violent relaxation discussed in Sec. 13.31 where we have to 
"guess" the form of C(f) that is relevant to our system (or try different functionals) . 

If we substitute Eqs. (I75 |l - (j73jl in Eq. (J7T|) and carry out the integrations on r" and v", we 
obtain 

df + Lf = e 3 r el^ f ds [ dVdVF'V -> r) t F"(r f -> r) t _ s 

1 d] 1 df' 







dt ' ~ J ~ r v dv" 

(75) 



C"(f) 9v u C"(f ) 9v' v 



We have written f t _ s = f(r'(t - s),v'(t - 8 ),t-s), f t _ s = f(r(t - s), v(t - s), t - s), F"(r' -> 
r) t = F»(r'(t) -> r(t)) and F u (r' -> r) t ^ s = F v (r'(t - s) -> r(t - s)) where r(t - s) and v(t - s) 
are the position and velocity at time t — s of the stellar fluid particle located in r = r(t), 
v = v(i) at time t. They are determined by the characteristics (J61j) of the smoothed-out 
Lagrangian flow. 

Equation (|75|) is a non Markovian integrodifferential equation: the value of / in r, v at time 
t depends on the value of the whole field f(r f , v', t — s) at earlier times. If the decorrelation 
time t is short, we can make a Markov approximation and replace the bracket at time t — s by 
its value taken at time t. Noting furthermore that the integral is dominated by the contribution 
of field stars close to the star under consideration (i.e. when r' — > r), we shall make a local 
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approximation and replace C"(f ) and J^r by their values taken at r. In that case, the foregoing 
equation simplifies in 



df -s _ 3 3 d 



ds I dVdViW -> rVF^fr' 



(76) 



x 



df 



C"(f')dv v C"(f)dv' v ) t 



where, now, / = f(r,v',t). The explicit reference to the past evolution of the system is only 
retained in the memory function 



J ds J d 2 r'F»(r' -> r) t F u (r' -> r 



t-s- 



This function can be calculated explicitly if we assume that, between t — s and t, the stars 
follow linear trajectories, so that v(t — s) = v and r(t — s) = r — vs [1. This leads to the 
generalized Landau equation 



(77) 

where K^ v is the tensor 
(78) 



dt dv 1 * 



f d 3 VK^\ — \j 
J \C"(f 



df 



07 



(f')dv u C"(f')9v' u j' 



= 2ttGW \nA- ( 5^ - 

u V u 



and u = v' — v, In A = ln(i2/e r ). This equation applies to inhomogeneous systems but, as a 
result of the local approximation, the effect of inhomogeneity is only retained in the advective 
term. 

By using a different approach based on a Maximum Entropy Production Principle (MEPP), 
Chavanis, Sommeria & Robert jH] have proposed a relaxation equation for p(r, v, rj, t) of the 
form 



(79) 



dp d r 



D 



dp 
dv 



From this equation, we can deduce a hierarchy of equations for the moments f n = j pi] n d 3 rd 3 vdr]. 
The equation for the first moment / is 



30) 



d -l + Ll-^{ 
dt +Lf ~ dv\ 



D 



<7V 



If we close the hierarchy of equations with Eq. ([74)1 . we obtain a self-consistent equation of the 
form 



(81) 



df +Ll- 9 { 



D 



df P ; 
dv C"(/) V 



Note that this equation can also be derived from Eq. (J77J) by replacing /' by its equilibrium 
value. Then, the diffusion coefficient D can be explicitly evaluated (see j2Sj in the two-levels 
approximation). Furthermore, Chavanis et al. [8 j proposed to let j3 depend on time, i.e. 
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(3 — so as to conserve energy. This heuristic "microcanonical description" is more adapted 
to the context of the violent relaxation than a "canonical description" with fixed /3. 

By developing a quasi- linear theory of 2D turbulence for the Euler-Poisson system [24], we 
have derived a kinetic equation of the form 



(82) ^ + uVu 

where K^ v is the tensor 
(83) 



/ ^ K ^m, a 



_/n 9u) du' 



K'"' — I <V"- 



and £ = r' — r. This equation is valid in the two-levels approximation of the statistical theory 
L) = {(To, 0}. Once again, the multi- levels case would lead to N coupled differential equations for 
p(r, a, t), the density probability of finding the vorticity level u = a in r at time t. Alternatively, 
if we implement a closure approximation of the form u 2 = l/C"(uf), similar to Eq. (|74p. we 
obtain a simpler equation 



54) 



du 
~dt 



+ uVu 



eV_9_ 
8n 2 dr^ 



d 2 r'K^ 



du 



1 du' 



C"(u') dr» C"(u) dr'' 



which takes into account several features of the process of "inviscid relaxation" in 2D turbulence. 
Unfortunately, as discussed in [21] , Eqs. (|82jl and ()84j1 do not conserve energy. Non-markovian 
effects may be necessary to restore the conservation of energy. 

By using a Maximum Entropy Production Principe (MEPP), Robert & Sommeria [5] have 
proposed a relaxation equation for the local distribution of vorticity p(r, er, t) of the form 



(85) 



g+nv, = v<i> 



From this equation, we can deduce a hierarchy of equations for the moments u r 
For the first moment uJ, we have 



J pa n da. 



(86) 



du 
~dt 



+ uVu = V<D 



VaJ + (3u 2 Vip 



Kazantsev et al. jzo] proposed to close the hierarchy of equations with a Gaussian approxima- 
tion. The resulting equation converges to a state of minimum enstrophy. More generally, we 
propose to close the hierarchy of equations by the relation u 2 = 1/C"(u) so that 



(87) 



du 
~dt 



+ uVu = V< D 



VaJ + 







C"(uJ) 



The function C(u) is a free function which has to be adapted to the context (see Sec. 13.3)1 . 
This indetermination reflects the fundamental observation that there is no universal form of 
entropy S\uJ] in 2D hydrodynamics [H]- The enstrophy |2E], the Fermi-Dirac entropy |2Zj, the 
Boltzmann entropy [2E] (leading to a sinh-Poisson equation) or the functional proposed by 
Ellis et al. ;6] (leading to a Gamma law for the vorticity fluctuations) are particular functionals 
S[u] which prove to be more relevant than others in a specific context (geophysical flows [2*3] . 
Jupiter's Great Red Spot 2D turbulence [301 and jovian atmosphere |E] respectively). None 
of them has a universal domain of validity but generalized entropies can possibly be regrouped 
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in "classes of equivalence" [3]. Tsallis entropies, sometimes arising in 2D hydrodynamics )3T| . 
are special because they are due to incomplete relaxation. They are particular H-functions 
associated to a dynamical equilibrium, not to a thermodynamical equilibrium (they are not 
true entropies). The self-confinement of 2D vortices can be explained by a lack of mixing |3*2*] . 
It can be taken into account in the statistical theory by using relaxation equations with a 
space-dependant diffusion coefficient related to the local fluctuations of the vorticity [33]. The 
same arguments can be invoked to account for the confinement of galaxies in astrophysics jH] . 

In conclusion, the generalized kinetic equations introduced previously can provide a simple 
parametrization of turbulence (mixing) in the context of the violent relaxation of stellar systems 
and 2D vortices. They can take into account incomplete relaxation thanks to a varying diffusion 
coefficient |3SIIS|- In addition, due to the thermodynamical analogy discussed in Sec. 13.2) they 
can also serve as powerful numerical algorithms to compute nonlinearly dynamically stable 
solutions of the Vlasov-Poisson or 2D Euler-Poisson systems. This is a great practical interest 
of these equations independently of the statistical theory. 



4 Properties of the generalized Landau equation 

4.1 Conservation laws 

In this section, we derive the conservation laws satisfied by the generalized Landau equation 
()2(J|) . Let us first introduce the current of diffusion 

= - / ^i^{*(/iW)|£ - s(/Wi)f|}- 

The time variation of energy can be written 

(89) E = f + rf3rrf3v = / J/ ' V rf3rrf3v ' 

where we have used an integration by parts. Introducing the current (J88|) in Eq. (J89|) . we 
obtain 

(90) E = - J dWvd^K^v^gifJhif)^ _ g ( f)h ( h )^l\ 
Interchanging the dummy variables v and vi, we get 

(91) E = J d\d\d*v x K^v^g{h)h{f)^ _ 

where we have exploited the symmetrical form of the diffusion current. Taking the half sum of 
the last two expressions, we find 

(92) E= l -J dWy^K-^J^H/)! - 9(f)Kfi)^y 
Noting that 

(93) K^u v = 0, 

according to Eq. 1)21)1. we finally establish that E = 0. Therefore, the generalized Landau 
equation conserves energy. We can show by a similar procedure that it also conserves angular 
momentum and impulse. 
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4.2 Generalized H-theorem 



The rate of production of generalized entropy is given by 

(94) S = - [ C\f)%d?vd\ = - [ C"(f)J»^-d 3 rd\ = - [ ^lr^Ldhd\, 
y ! J y 1 dt J KJ) dv» J g(f) dv» 

where we have used Eq. to get the last equality Inserting the current (fKKJ) in Eq. (jMJ), 
we obtain 

Interchanging the dummy variables v and v 1; we get 

Taking the half sum of the foregoing equations we find that 

(97) ^l/iWv^^-^}^^-^}. 

Noting that X^K^X" = X 2 - (X ■ uf/u 2 > 0, we conclude that S > 0. Therefore, the 
generalized Landau equation satisfies a generalized if -theorem. 

4.3 Equilibrium distribution 

Taking the derivative of Eq. (|T7j) with respect to v, we get 

(98) c"(/ e? )^ = -/?(v-nxr-U). 
Inserting this relation in the current (J88|) and using Eq. (J29|) . we obtain 

(99) Jf = -(3 J d\ 1 K^gg l {v\ - v"), 

which vanishes identically in virtue of Eq. (|93p. Therefore, an extremum of the generalized 
entropy at fixed mass, energy, angular momentum and impulse is a stationary solution of the 
generalized Landau equation. Alternatively, a stationary solution satisfies S = hence J/ = 0. 
We can show that this condition implies that f eq is of the form ()27)) . The proof is similar to 
the one given for the ordinary Landau equation. 

Finally, considering the linear stability of a stationary solution of the generalized Landau 
equation, we can derive the general relation |3] 

(100) 2X5 2 J = 5 2 S > 0, 

connecting the growth rate A of the perturbation (such that 5f ~ e xt ) to the second order 
variations of the free energy J = S — (3E and the second order variations of the rate of entropy 
production 5 2 S > 0. This aesthetic relation implies that a stationary solution of the generalized 
Landau equation is linearly stable if, and only if, it is an entropy maximum at fixed mass and 
energy [3]. 
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5 The thermal bath approximation 



The generalized Landau equation (J2Uj) can be put in a form reminiscent of a Fokker-Planck 
equation 



(101) 



df d 



dt dv" 
by introducing a diffusion tensor 



(102) 

and a friction term 
(103) 



V 



= I K^g{f x )d^ 



The ordinary Fokker-Planck equation [2] can be written 

df 



(104) 
where 
(105) 



d 2 d 
_ f /r/^ + —(ffi*) 



c 



1 (At> ^Af v ) 

2 At 



(Aw") 

~aT' 



are the first (friction) and second (diffusion) moments of the velocity deviation. An equivalent 
form of the Fokker-Planck equation is 



(106) 



d L = A_ 

dt dv» 



dv u 



For the ordinary Landau equation, we have g(f) = f et h(f) = 1. Therefore, Eq. (jlUlj) 
becomes 



(107) 



df 







df 



dt dv» 

By comparing with the Fokker-Planck equation (|l()fi|) . we obtain the classical relations 



(108) 



+ C M = rf. 



According to Eq. ()108j) . there is a difference between ff and (Av) due to the velocity dependance 
of the diffusion tensor D^ v ' . Noting that 



(109) 



dv u J dv 
we find that Eq. (ITU5J) yields 

(110) 



fid\x 



l(Av) 



2 At 

Therefore, in the classical case, the vector — rj represents half the friction force. 



20 



The generalized Landau equation (jl(Jl|) can be rewritten 

where if is a primitive of h. By comparing with the ordinary Fokker-Planck equation ()104|) . 
we get 

If we impose the relations g(f) = f and h(f) = fC"(f) of Sec. \'2A\ Eq. (jll2|) can be rewritten 

c(f)' 



1 (Av^Av u ) 

113 - v ' = D^f 

v ' 2 At J 



f 



This relation was noted in In the particular case of the Tsallis entropy, Eq. (|113|) reduces 
to the form considered by Borland In the Borland approach, generalized Fokker-Planck 
equations arise because the transition probabilities ^ u and C M depend explicitly on the distri- 
bution function. A different approach is followed by Kaniadakis JU]- Starting from a kinetic 
interaction principle (KIP), he obtains a generalized Fokker-Planck equation of the form 



(114) 9 



dt dv» 



This corresponds to his Eq. (17) in our notations. Comparing with our equation (jlOlj) . we 
find that the classical relations (|108|) are preserved in this generalized framework. This is much 
more aesthetic than Eq. ()112|) . This is also more physical because we would expect that the 
transition moments depend only on the form of interaction, not explicitly on the distribution 
function. Besides, the other approach leads to ambiguity because, comparing Eqs. (jl(Jl|) and 
(HHD, we get 

(115) C u = D^h(f), Ji_ + C m=^V, 

which differ from Eq. (J 112)) except in the classical case. 

Physically, the Landau equation describes the statistical evolution of a system of self- 
interacting particles in a mean-field approximation. This description treats all the particles 
on the same footing and conserves the energy of the whole system. This corresponds to a 
microcanonical description. We can alteratively view the Landau equation as describing the 
evolution of a test particle (described by /(v,t)) interacting with field particles (described 
by /(vi,t)). In the so-called thermal bath approximation (canonical description), we consider 
that the field particles are in statistical equilibrium and replace /(vi,t) by their equilibrium 
distribution / eg (vi) given by Eq. (|27|). This approximation transforms an integro-differential 
equation (Landau) into a differential equation (Kramers). Combining Eqs. (|1(J3|) . (|98|) and 
(EH, we get 

(116) rf = pj ^(AKd 3 ^, 

with w = v — fixr — U. Now, using Eqs. and (jl(J2|) . we have equivalently 

(117) rf = [ K^gi^w'd 3 ^ = (3D^w v . 
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This relation can be interpreted as a generalized Einstein relation. In the test particle approach, 
the kinetic equation (|101|) thus becomes 



;ii8) 



dt 



d 



D 



This equation will be called the generalized Kramers equation. Assuming for simplicity that 
the diffusion is isotropic, taking Q = U = and imposing the relations of Eq. ({HI} , we recover 
the generalized Kramers equation 



(119) 



dt dv\ 



proposed in [Hj. In the present approach, the generalized Kramers equation is derived from 
the generalized Landau equation by a systematic procedure. This makes possible to determine 
explicitly the expression of the diffusion coefficient (see Sees. |3|7J) which was left unspecified 
by the MEPP approach 0. 

We recall also that in the long time limit (or strong friction limit £ = D/3 — > +oo), the 
generalized Kramers equation leads to the generalized Smoluchowski equation 

dp 



(120) 



dt 



(Vp + pV$) 



where p(p) has the interpretation of a pressure [3 . The generalized Smoluchowski equation was 
first introduced in jSJ. The gravitational (i.e. attracting) Smoluchowski-Poisson and generalized 
Smoluchowski-Poisson systems have been studied in [331 EH! (see also jSZj for a connexion 
with mathematical results). We shall say that a kinetic equation has a generalized canonical 
thermodynamical structure if it increases the generalized free energy J = S—f3E at fixed inverse 
temperature (3 and mass M (the function J differs from the usual free energy F = E — TS 
by a factor —(3). We have shown that the generalized Kramers equation and the generalized 
Smoluchowski equation possess such a thermodynamical structure |3. 3 . We can formally obtain a 
microcanonical formulation of the generalized Kramers and Smoluchowski equations by letting 
the inverse temperature (3 depend on time, i.e. (3 = (3(t), so as to conserve energy. With 
this modification, the generalized Kramers and Smoluchowski equations have a microcanonical 
thermodynamical structure 



6 Diffusion coefficient and dynamical friction 
6.1 Generalized Rosenbluth potentials 

In order to evaluate the expressions of diffusion coefficient and friction, it is convenient to 
introduce auxiliary functions which are called the Rosenbluth potentials Noting that 



121) 



d 2 u 
dv^dv" 



u 



2 5^ - uV 



we can rewrite the diffusion tensor (jl()2j) in the form 



(122) 

where 

(123) 



= A 



d 2 X 



(v), 



X(v) 



g(h)W- vi|d 3 vi, 
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is the generalized Rosenbluth potential associated to the diffusion. Writing the friction term 
(J103J) in the form 

(124) = - J K^f^Mtfy^ 



and integrating by parts, we obtain 

= ('. 

dv u x 



(125) rf = I ff(A)dV 



Noting that 

QKvv oi u (9/1 

(126) = 2 A = -2A 1 



we can rewrite the friction term as 
(127) rf = -2A—M, 



where 



(128) A(v) = / -^L^-d^, 

J | v — v i| 

is the generalized Rosenbluth potential associated to the friction (we need to impose that 
H(fi) tends to zero for |vi| — > +oo to make the integral well defined; this fixes the constant of 
integration in H). We note also that 

(129) — — = / - in -g(fi)d 3 v 1 = - / ^-^(/i)d 3 vi = 2A- 



dv v J dv u J dv\ dv^ 

where 

(130) a(v) = / d 3 Vl , 

J | v _ v i| 

is the potential associated to the velocity dependance of the diffusion coefficient. In the usual 
thermodynamical framework where g(f) = H(f) = f, the potentials A and o coincide. 



6.2 Isotropic distribution of velocities 

When the velocity distribution of the field particles is isotropic, i.e. j\ = /(i>i), we can obtain 
more explicit expressions for the Rosenbluth potentials. For reasons of symmetry, the Rosen- 
bluth potentials depend only on v — |v|. To determine x( v ) an d A(t>), we use the identity 



(131) 




where t>< and v > denote the smallest and largest value of v and v±, Pi{x) is a Legendre polyno- 
mial and 7 is the angle between v and Vi- Then, we have 

+~ r+oo 2 1 P7T 

(132) X(v) = 2ny2 -^H[fi^)\dvi ^(cos 7 ) sin 7 rf 7 . 
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With the change of variable x = cos 7, we get 

+00 „ +00 2.A 



(133) X(v) = 2tt V / J^Hlf^dv! / Pi{x)dx. 

1=0 J ° v > J- 1 

Using the identity Jj^ 1 Pi(x)dx = 2<5 ro , the foregoing expression can be simplified in 

p+00 2 |"i rv r+oo 

(134) \(v) = 4tt / -^Hif^dvx = 4vr - / v'ff [/(vx)]^ + / v x H[f{v x ))d 

JO v > l v JO Jv 



Similarly, we have 

(135) a(v) = 4tt 



vig[f {vi)}dvi 



To determine x{ v )i we write |v — Vi| = {v 2 +v \ — 2vv\ cos7)/|v — v x | and we use the identity 
ffTTTTj). We then obtain 



(136) 



+00 „+oo 2 / 



i=0 



/*7T 

fl , [/(' u i)]^ l 'i / ( t ' 2 + v i ~ %vv 1 cos7)Pi(cos7) sin 7^7. 
Jo 



(v 2 + v\ — 2vv\x)Pi(x)dx. 



With the change of variables x = cos 7, we get 

+00 „ +00 21 r+i 

(137) X (v) = 2irJ2 ~1^9[f{^)\dvi ( 

;=o Jo v > J- 1 

Using the identity Jj^ 1 xPi(x)dx = we obtain after some manipulations 

(138) X (v) = — 



6.3 The diffusion coefficient 

We are now in a position to determine an explicit expression for the diffusion tensor D^ u that 
is valid for an arbitrary isotropic distribution function of the field particles. Starting from the 
identity 



(139) 



d 2 x v»v u (d 2 x ldx\ , Idx,^ 



dv fJj dv u v 2 \ dv 2 v dv J v dv 
the diffusion coefficient (|122j) can be put in the form 



(140) 

where 
(141) 
and 
(142) 



D 



D „ - ~ 2 D, 



v^v u 1 

^4- + -D ± 8»\ 



v dv 



Du = A 



d\ 

dv 2 ' 



24 



are the diffusion coefficients in the directions perpendicular and parallel to the velocity of the 
test particle. To see that, we consider a system of coordinates where the z-axis is taken in 

v. In this system of coordinates, all the 



the direction of v so that v x = v y = and v 
non-diagonal elements of D^ u vanish while D X1 
formula (|138|) . we have explicitly 



D 



yy 



\D< L and D z 



Du. According to 



(143) 



8tt „ 1 
— A- 

3 v 



I \ 3vl-^g[f(v l )]dv 1 + 2v jf vig[f(vi)]dvi 



(144) 



87r ,1 

D \\ = T A ~ 
3 v 



" v 4 

-^g[f{vi)]dvi + v 
v z 



vig[f(vi)}dv 1 



6.4 The dynamical friction 

The friction term ()127|) can be simplified similarly. For an isotropic velocity distribution, 

d\ dX 



(145) 
so that 
(146) 

Using Eq. (|134|) . we get 
(147) 

Similarly, we have 
(148) 



dv^ dv v 



ri = -2A — — v. 



v dv 



r) =87tA^ f vfHlfiv^dvx. 
v Jo 



Q Dii V 

dv v 



-8vrA^ f tf x g[f{v x )\dv x . 
v Jo 



We note that when the field particles have an isotropic velocity distribution, the dynamical 
friction experienced by the test particle 



(149) 



(F) 



friction 



-8nA— I v, 
V s Jo 



f 



f 



dvi, 



is parallel and opposite to its velocity v. Moreover, it is due only to field particles with velocity 
v\ < v. In fact, we can obtain this result without calculation. Indeed, according to Eqs. (|127j) 
et (ED, we have 



(150) 



■q = 2A 



V — V X 
|v — Vi|' 



H{fi}d 3 vi, 



and a similar expression for dD^ v / dv v . Now, Eq. ()150j) is analogous to the gravitational force 
created in v by a distribution of mass with density H[f(vi)} where v plays the role of the 
position r. According to Newton's law, the gravitational force created in r by an isotropic 
distribution of mass depends only on the mass interior to r and is given by an expression 
equivalent to Eq. ()147|) . 
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7 Explicit results for typical distribution functions 



7.1 Isothermal distributions 

We shall now obtain explicit expressions of diffusion coefficient and friction force for typical 
distribution functions. First, we consider the case of ordinary thermodynamics based on the 
Boltzmann entropy. In the thermal bath approximation, the field stars are at statistical equi- 
librium with the isothermal distribution 



(151) 



/e,(v 1 )=p(^) 3/2 e-4. 



For the Boltzmann entropy C(f) = / In/, we have g(f) = H(f) = f . The diffusion coefficient 
and the friction coefficient can then be calculated explicitly by substituting Eq. (|151j) in Eqs. 
(|143j) . (|144j) and (|147|) . and performing the integrals. Introducing the notation X = j3v 2 /2 
and the function 



(152) 

where 

(153) 



G(X) 



2 1 



X 2 



x 



re~* dt 



2X 2 



x 



erf(X) 



2X 



erf(X) = -= 



"* dt 



is the error function, we find after some elementary calculations that 



(154) 



2ApG(X)~, 
v 



(155) 



D± = 2Ap[erf(X) - G(X)}-, 

v 



(156) 



y=2Ap/3G(X)~. 



These results are well-known and are recalled here for sake of completeness. We note that 
D\\ ~ v~ 3 for |v| — > +oo. On the other hand, combining Eqs. (jllOj) . (jl56j) and p54|) we note 
that the dynamical friction can be written in the form 

(157) (F) friction = -2n = -2/3D\\v, 

so that the Einstein relation reads £ = 2/3 D\\ (the factor 2 is due to the velocity dependance of 
the diffusion coefficient). Quite generally, in the thermal bath approximation, we have 



(158) ri 
which results from Eqs. flll7j) and (jl4()|) . 
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7.2 Fermi-Dirac distributions 



We now consider the case of quantum particles (fermions) described by the Fermi-Dirac entropy. 
In the thermal bath approximation, the field particles are at statistical equilibrium with the 
Fermi-Dirac distribution 



(159) 



1 + Ae^ 

The parameter A is related to the density p by 



(160) 



where 



(161) 



P 



In(t) 



£3/2 



-WA), 



+oo 



X 



l + te 1 



-dx, 



is the Fermi integral of order n. For the Fermi-Dirac entropy, the functions g and H are given 
by g(f) = fi^ — f/Vo) an d H(f) = f. The diffusion coefficient and the friction coefficient can be 
calculated explicitly by substituting Eq. ()159|) in Eqs. ()143j) . ()144j) and ()147j) . and performing 
the integrals. Introducing the notation X = a/ '(3v 2 /2 and the incomplete Fermi integral 



(162) 



In(t,X) 



hit) 



X 



X 



1 +te x 



-dx, 



we find after elementary calculations that 



(163) 



D\\=Ap— i 1/2 (x,xy 



(164) 



D ± = Ap 



^^/_ 1/2 (A,X)-i I / 1/2 (A,X) 



r, = Ap^I l/2 {\X)- = (3D\\v 
X z v 



(165) 



These equations were previously derived in [zo] in a less neat form. For t — > +oo, we have the 
equivalent 



(166) 

On the other hand, 
(167) 

Therefore, 
(168) 



I n (t) ~ ir(n + i). 



J n (+oo,X)= 2 / y 2n+1 e-y 2 dy. 
r(n + 1) Jo 



h /2 (+oo, X) = 2X 2 G(X), /_i/ 2 (+oo, X) = erf (X) 



In the non-degenerate limit A — > +oo, Eqs. ()163j) - (jl65j) return the results p54|) ~ p56j) valid for 
classical particles. 
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7.3 Tsallis distributions 



We finally consider the case of generalized thermodynamics based on Tsallis entropy. In the 
thermal bath approximation, the field stars are at statistical equilibrium with the polytropic 
distribution 



(169) 



f eq (y l ) = A[\-j- 



if v < v m = v2A and / = otherwise (we restrict ourselves to q > 1). The polytropic index n 
is related to the g-parameter [7j by 

3 1 
2 + q-1 



(170) 



n 



The parameters A and A are related to the generalized temperature f3 and to the density p by 
the relations 



(171) 

where 

(172) 



.4 



q-l 



P = 4nV2A\ n B( l,n- I 



B(m, n) 



_m— 1 



(1 - x) n - l dx, 



is the beta-function. For the Tsallis entropy, the g and H functions are g(f) = f and H(f) = f q . 
The diffusion coefficient and the friction coefficient can be calculated explicitly by substituting 
Eq. ()lfi9j) in Eqs. (|14H|1 . (j!44|) and p47j) and performing the integrals. Introducing the notation 
X = \J{n + l)t> 2 /2A and the incomplete beta-function 



(173) 



B x (m,n) 



B(m,n) Jo 



71+1 



if X < \Jn + 1 and Ix( m , n) 
(174) 



1 otherwise, we find after elementary calculations that 



(175) 



D ± = Ap 



(176) 

For n — > +oo, we have the equivalent 



,= ^,(J, n + i)^l,v. 



(177) 

On the other hand, 
(178) 

Therefore, 
(179) 



B(m, n) 



T(m) 



B x (m, +oo) 



T(m) 



y 2m-l e -y- dy 



fix [ ~ +oo j = erf (X), B x f^, +oo ) = 2 A' 2 6'( X). 



In the limit n — > +oo, corresponding to g — > 1, Eqs. ()174j) - (jl76j) return the results f)154[) - f|156^ 
obtained in the context of ordinary thermodynamics. 
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8 Truncated models accounting for a permanent escape 
of particles 

In this section, we shall derive generalized truncated distribution functions accounting for an 
escape of particles above a limit energy e m . For simplicity we shall work in the thermal bath 
approximation. We thus describe the evolution of the distribution function by the generalized 
Kramers equation 



180) 



dt dr 



<9v 



d_ 
<9v 



K 

1 ! 3 



h(f)^ + Pg(f> 



where we have neglected anisotropy, taken D^ v = D^S^ in order to respect the Einstein relation 
(j!58|) and used the asymptotic expression Du ~ v~~ 3 of the diffusion coefficient valid for high 
velocities. We assume that high energy particles are removed by a tidal field (for globular 
clusters, this is the gravitational attraction of a nearby galaxy). We seek therefore a stationary 
solution of equation (jl80j) of the form / = /(e) satisfying the boundary condition f(e m ) = 0, 

2 

where e = y + $ is the energy of a particle and e m is the escape energy above which / = 0. 
Using the identity ^r(-r) = (valid for large |v|), we obtain 



(181) 

or, equivalently, 
(182) 



T h(f)f + /3g(f) =0. 
ae ae 



h(f)f + Pg(f) = -J, 
ae 



where J is a constant of integration representing physically a current of diffusion. If J = 0, the 
foregoing equation reduces to 



183) 



C"(/)| + /3 = 0, 



where we have used Eq. ()29|) . After integration, we recover the usual equilibrium distribution 
(184) c'(f eq ) = -(3e-a. 



If J 0, Eq. ()182j) accounts for an escape of particles at a constant rate J. The system is 
therefore not truly static since it looses gradually particles but we can consider that it passes 
by a succession of quasi-stationary states that are solution of Eq. ()182|) . Equation ()182|) is a 
first order differential equation which can be integrated as 



185) 



h{t)dt 
o g{t) + J//3 



/3(e m - e). 



Let us consider particular cases of this equation. For classical particles described by the 
Boltzmann entropy, Eq. ()185|) reduces to 



186) 



/ dt 



Pi^m ~ e). 



/o t + J/0 

The integral is readily performed and we obtain the Michie-King model 
(187) / = A{e~ Pe -e~ Ptm ), 
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where we have set A = (J/ft)e l3em . The Michie-King model describes the tidally truncated 
structure of globular clusters in astrophysics [TT]. For quantum particles described by the 
Fermi-Dirac entropy, Eq. (j!82|) reduces to 



(188) 



dj_ 

de 



+ Pf(i-f/m) = -J- 



This is a Riccatti equation that can be solved analytically. Assuming that degeneracy is negli- 
gible for energies close to the escape energy, we get 



(189) 



/ = Vo- 



-Per, 



This truncated distribution was previously derived in [23]. It could describe the case of galactic 
halos (e.g., massive neutrinos in dark matter models) limited in extension by tidal forces. It 
could also be of interest for collisionless stellar systems with Lynden-Bell's interpretation of 
degeneracy. Finally, for systems of particles described by the Tsallis entropy, Eq. ()185|) reduces 
to 



(190) 

For q = 1, we recover Eq. 
(191) 

we can rewrite Eq. 
(192) 



W J V- X dx {J 



l-q 



ft (em - e). 



1 + * q 
. If we introduce the function 

X t q-l 



<j> q (x) 



l + t 



-dx, 



in the form 
J 



./ 



ft 



(J/ft) 



l-q 



Q 



ft(e r 



For q — 1, we recover the Michie-King model. A detailed study of these truncated models will 
be given elsewhere. 



9 Conclusion 

In this paper, we have shown that standard kinetic equations (Boltzmann, Landau, Kramers, 
Smoluchowski,...) can be generalized so that they satisfy a i?-theorem for an arbitrary func- 
tional of the form S = — J C(f)d 3 rd 3 v, where C(f) is a convex function. Boltzmann, Fermi, 
Bose and Tsallis entropies are particular functionals of the above form. These generalized 
kinetic equations have a thermodynamical structure which corresponds either to a microcanon- 
ical (Boltzmann, Landau) or a canonical (Kramers, Smoluchowski) description. One important 
conclusion of our work is that Tsallis entropy does not play any special role in this generalized 
thermodynamical formalism except that leading to simple distributions [3]. Therefore, the ques- 
tion that naturally emerges is whether other arguments can give Tsallis entropy a fundamental 
justification or whether Tsallis entropy is just a simple functional (associated to power-laws) 
extending Boltzmann entropy and providing a good fit of several observed phenomena. This is 
clearly an important point to be settled in the future. 

There are at least two distinct notions of generalized thermodynamics. Generalized thermo- 
dynamics and kinetic equations can arise in complex media when the transition probabilities 
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have an expression different from the one we would naively expect |lL)j . We believe that gen- 
eralized kinetic equations are essentially effective equations attempting to take into account 
"hidden constraints" that are not directly accessible to the observer or that are difficult to for- 
malize. If this idea is correct, it means that it will never be possible to justify these equations 
from first principle (except in toy models). This also explains naturally why there is some 
indetermination in the theory, either in Tsallis g-parameter or more generally in the function 
C(f). The important point in that context is to understand why standard kinetic theory breaks 
up. To our point of view, it is not sufficient to say "since the system is non-extensive, Boltz- 
mann entropy is not correct and Tsallis entropy must be used instead" . This argument is too 
simplistic and does not bring any physical insight in the problem. The main interest of Tsallis 
nonextensive thermodynamics (in addition to the development of a consistent formalism) is to 
show that the naive Boltzmann description fails in many systems. However, the reason of this 
failure has to be understood in each case and this is the main challenge. 

A notion of generalized thermodynamics also emerges in the context of the violent relaxation 
of stellar systems and 2D vortices (or other Hamiltonian systems with long-range interactions) 
described by Vlasov-type equations. In that context, we can explicitly illustrate the notion of 
"hidden constraints" that we mentioned previously. According to rigorous statistical mechan- 
ics |18| I19j. the metaequilibrium state resulting from complete violent relaxation is obtained 
by maximizing the Boltzmann entropy S[p] for the fine-grained distribution p(r,v,i]) while 
conserving mass, energy and an infinity of additional constraints played by the Casimirs. It 
turns out that the coarse-grained distribution function /(r,v) also maximizes a functional 
S[f] — — J C(f)d 3 rd 3 v while conserving only mass and energy (robust constraints) [HIS]- This 
functional, which could be called a generalized entropy, is non-universal due to fine-grained 
constraints that depend on the initial conditions. The Casimirs represent "hidden constraints" 
because, in practice, we just know the coarse-grained field and we do not have access to the 
initial conditions. In case of complete violent relaxation S[f] is never Tsallis entropy since 
/ > according to Lynden-Bell's theory. However, violent relaxation is incomplete in general. 
The only thing we know for sure is that the metaequilibrium state reached by the system is a 
nonlinearly dynamically stable stationary solution of the Vlasov equation on a coarse-grained 
scale. This implies in many cases that it maximizes a if-function S[f] — — J C(f)d 3 rd 3 v at 
fixed mass and energy [7]. There, C(f) depends on the initial conditions and on the efficiency 
of mixing. Tsallis entropy is a particular H-function. Since the dynamical stability criterion is 
similar to a generalized thermodynamical criterion, we believe that it is this thermodynamical 
analogy that justifies the consideration of Tsallis functionals (they are not true entropies!) in 2D 
turbulence and stellar dynamics. However, we stress that Tsallis entropy has no fundamental 
justification in that context [32] and that, indeed, most vortices and galaxies are not polytropic 
012]. Finally, we emphasize that the true statistical equilibrium state resulting from encoun- 
ters between stars or between point vortices (collisional relaxation) is described, in a suitable 
thermodynamic limit, by the ordinary Boltzmann entropy although the system is non-extensive 
and non additive. The pecularities due to the absence of entropy maximum in gravitational 
systems correspond to important physical processes (evaporation and gravothermal catastro- 
phe) and not to a break up of thermodynamics [7]. We hope that this critical discussion will 
help to clarify the different notions of "generalized thermodynamis" that appeared in the recent 
literature. 

I acknowledge stimulating discussions with F. Bouchet, T. Dauxois and C. Sire. 
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